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We show that the application of twisted mass QCD (tmQCD) with four (Wilson) quark flavours to the compu- 
tation of lattice weak matrix elements relevant to AI ~ 1/2 transitions has important advantages: the renormali- 
sation of A' — > tt matrix elements does not require the subtraction of other dimension six operators, the divergence 
arising from the subtraction of lower dimensional operators is softened by one power of the lattice spacing and 
quenched simulations do not suffer from exceptional configurations at small pion mass. This last feature is also 
retained in the tmQCD computation of A" ^ tttt matrix elements, which, as far as renormalisation and power 
subtractions are concerned, has properties analogous to the standard Wilson case. 



1. Introduction 

At scales well below Mw, but above the charm 
quark mass, the effective weak Hamiltonian for 
CP-conserving, AS* — 1 decays can be written 
as: 

He// - K<iK;^[c+ WMw)o+(m) 

+ C_WMh')0«(m)] , (1) 

where 0^(/x) are the dimension-6 four-fermion 
operators 

— Oyy^j^^ — Oyj^^j^y (2) 

renormalised at a scale /i and 7^ = 57;i(l — 75) 
(the subscript R indicates renormalised quanti- 
ties). Parity ensures that — > tt+tt" and 

—> matrix elements will receive contribu- 
tions only from Oyj^^j^y, while tt'^ matrix 
elements will arise only from Oyy^j^j^. 

Under isospin transformations, the operator 

is purely / = 1/2, while exhibits both 
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/ = 1/2 and / = 3/2 parts. The A/ = 1/2 rule 
can be accounted for by an enhancement of the 
K — > TTTT matrix element of C-0~ with respect 
to that of The contribution to this en- 

hancement coming from the ratio of Wilson coef- 
ficients C_/C+ is small in typical renormalisation 
schemes (such as MS); thus, the bulk of the en- 
hancement must come from the ratio of the ma- 
trix elements of the operators themselves. 

2. Wilson fermion renormalisation 

A serious difficulty arises upon attempting to 
extract physical amplitudes directly from Eu- 
clidean correlation functions corresponding to 
matrix elements with more than one particle in 
the final state [Q. Although it has recently been 
shown that the problem can be bypassed , there 
is still scope in adopting the long-standing alter- 
native to a direct calculation. This consists in us- 
ing chiral perturbation theory to obtain K ^ tttt 
matrix elements from if — > tt ones the latter 
not being plagued by the presence of final state 
interactions. 

We are then faced with the problem of oper- 
ator renormalisation and mixing. Besides a log- 
arithmically divergent renormalisation constant, 
the operators also mix with two dimension-3 
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operators {sj^d and sd with coefficients Cp and 
respectively) and two dimension-5 operators 
[sa^yGfi^d and sa^uG^j^ud with coefficients 
and respectively). In the case of Ginsparg- 
Wilson fermions, chiral symmetry ensures that 
these subtractions are mild, since Cpg ~ 
and df ^ ^ C'(a^) 0- Once chiral symmetry is 
lost with Wilson fermions, renormalisation pat- 
terns significantly worsen. Ignoring operators 
with contributions to matrix elements which van- 
ish by the equations of motion we now have (on 
the basis of CPS symmetry) 



+Cp{gl,m,a) s-fsd 
+df{glm,a)sa^^F'"'d^ 



^VA+AV 



, (3) 



'^VV+AA 



k=l 

+dtiglm,a)sa^,F''''d+---] , (4) 



where k ^ VV - AA,SS,PP,TT (in stan- 
dard notation) and the ellipses indicate subtrac- 
tions which, being 0{a), are unimportant to the 
present discussion^. Three problems immediately 
arise. First, the subtractions of the dimension- 
3 operators involve a linearly diverging coeffi- 
cient in the case of Ova+av and, even worse, 
a quadratically diverging coefficient for Ovv+aa'- 



- (mc - rriu) {rris ~ mi) 



1 



(5) 



A second, milder shortcoming is the mixing of 
the parity-even operator Oyy^^^ with four other 
dimension-6 operators (just as in the more famil- 
iar case of Bk)- Finally, Wilson fermions in the 
quenched approximation arc plagued by excep- 
tional configurations. In rcf. it has been shown 
that the tmQCD formulation of Wilson quarks 
does not suffer from the last two problems. Here 

^The coefficient in eq. (^) is also 0[a) and could have 
been omitted. 



we will demonstrate that the implementation of 
tmQCD can also reduce the quadratic divergence 
Cii K -K matrix elements to a linear one. 

3. tmQCD with four quark flavours 

For a recent review of lattice tmQCD see ref. § 
and references therein. Here we extend the for- 
mulation to four quark fiavours. The action is 
given by 

+ iphip + rrih + ifJ.h^5)il^h ■ (6) 

We distinguish a light quark doublet tpf — {u, d) 
and a heavy one ■0^ = {s,c). The 2x2 quark mass 
matrices are mi^h — diag(r7i„_s , md,c) and the 
twisted mass matrices are ni^h = diag(^u_s , /id,c)- 
For simplicity we impose mass degeneracy in the 
light sector; i.e. iriu = rud and Hu = —fJ-d- Two 
twist angles are then defined through ratios of 
renormalised mass parameters: 



tana 



tan/3 : 



(7) 



The equivalence of tmQCD and standard QCD 
is formally established through the axial field ro- 
tations 

^3- 



exp 



exp 



■0; : 
■iph 



(8) 



(and similarly for ipiji). 

For the iiT — > tt weak matrix element a conve- 
nient choice is given hy a = (3 = 7r/2. Operators 
are then related as follows 



[Sk 

['^VA+AV 



tmQCD 



I tmQCD 



I tmQCD 



QCD 



QCD 
QCD 



E US and Pk = U75S. 



(9) 



where = d-y^u, Sk = us and Pk = U75S. The 
following equation between tmQCD and standard 
QCD renormalised correlation functions (at un- 
equal space-time arguments) is thus true up to 
discretisation effects: 



{P^ix)O$^^^y{0)SKiy)) tmQCD = 
{P^ix)O^y^^^{0)PKiy)) QCD 



(10) 
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In the asymptotic limit the r.h.s. yields 
{n\Oyy_^_j^j^\K) QCB', this matrix element can thus 
also be computed in tmQCD (the l.h.s). We shall 
show below that in the latter formalism the renor- 
malisation properties of the four-fermion operator 
are much more convenient. 

For the K — > tttt matrix element a convenient 
choice of twist angles is a — —f3 = 7r/2. The 
equation between the corresponding renormalised 
correlation functions (at unequal space-time ar- 
guments) is then 

(P.(x)Pi(y)O±^+^^(0)5K(z))tmQCD 

(11) 

where now Pk = ^753, Sk = ds. In this 
case there are no advantages to be gained in the 
tmQCD computation as far as operator renormal- 
isation properties are concerned. However, in the 
quenched approximation, the tmQCD computa- 
tion is free of exceptional configurations. This is 
an important advantage close to the chiral limit. 

4. tmQCD renormalisation 

For the "convenient" choices of twist angles (i.e. 
a = 7r/2,/3 = ±7r/2) and using the discrete sym- 
metries of the tmQCD action (cf. we find 
the following renormalisation pattern for the op- 



erators O 



± 

VA+AV 



{ORjvA+AV - ZvA+AV 



'^VA+AV 



- sd - 



(12) 



where the leading behaviour of the coefficients (in 
the chiral expansion) is 



If 

a 



p{9o) (/"c 



^hi) (ps - fj-d) 



■Gsido) (Aic-/^«) (m^ 



(13) 



Thus, we now have a linear divergence in Cp 
(while contains two terms, both of 0(1)). This 
situation compares favourably to the standard 
QCD case, characterized by a quadratic diver- 
gence {cg ~ 1/a^; cf. eq. (||)). Moreover, there 



are no dimension-6 operators to be subtracted in 
the tmQCD case (in standard QCD there are four 
such subtractions; cf. eq. (^). 

In order to determine the linearly divergent co- 
efficient Cp we must resort to parity restoration 
in the continuum limit. Parity is broken by the 
lattice (Wilson) tmQCD but is recovered after 
renormalisation. In particular for the twist angle 
values of a = /3 = 7r/2, parity transformations in 
the continuum limit assume the form 

u{x) Z7o75 u{x) d(x) -i-inlb d(x) , 
s{x) «7o75 s{x) c{x) -17075 c(i) , (14) 

with X = (xq,— x) and similarly for the anti- 
quarks. Thus Oyj^j^j^y is a positive parity eigen- 
state whereas Pds = d'^^s is a negative parity 
eigenstate. This implies that for x ^ Q we have 

{[Or] VA+Avi^) iPR)ds ix))tmQCB — . (15) 

Expressing the above in terms of bare correlations 
yields 



C^VA+AV - 



+ C${PsdPds) 

+ 4{SsdPds)=0- (16) 

Analogous parity arguments can be used to show 
that the last term of the above expression is a 
lattice artifact (i.e. Cg (SgdPds) — C>{a)) and may 
be dropped. Thus, eq. ( p^ ) can be solved to de- 
termine Cp. 
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